Surface relaxation and ferromagnetism of Rh (001) 
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The significant discrepancy between first-principles calculations and experimental analyses for the relaxation of 
the (001) surface of rhodium has been a puzzle for some years. In this paper we present density functional 
theory calculations using the local-density approximation and the generalized gradient approximation of the 
exchange-correlation functional. We investigate the thermal expansion of the surface and the possibility of surface 
magnetism. The results throw light on several, hitherto overlooked, aspects of metal surfaces. We find, that, when 
the free energy is considered, density-functional theory provides results in good agreement with experiments. 
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The significant discrepancy between first-principles 



calculations ||1| 
(LEED) analyses 




and low-energy electron diffraction 
for the relaxation of the (001) 
surface of rhodium has been a puzzle for some years. The 
earlier LEED studies ||,0] concluded that the interlayer 
spacing of the surface layer {di2) is nearly identical to 
that in the bulk (do), he., the top-layer relaxation was 
determined to be Adi2/do = +0.5 ± 1.0 %. A recent 
LEED study § found Ad^/do = -1.16 ±1.6 %. On the 
other hand, first-principles calculations showed a large 
top-layer relaxation ranging from —3.2 % to —5.1 %, de- 
pending on the calculational scheme and/or the employed 
numerical accuracy [^J|,^,|| . Inward relaxations are in- 
deed the expected behavior of transition metals surfaces 
(see e.g. Ref. M ), and the practical zero relaxation de- 
termined by LEED is at least unexpected. 

In order to reconcile this disagreement between their 
calculations and experiment, Feibelman and Hamann |l|] 
proposed that in the experimental study the metal sur- 
face may be contaminated by residual hydrogen adsorp- 
tion (see also Ref. |^). Indeed, hydrogen is not easy 
to detect and quite soluble in transition metals, such as 
Ru, Rh, and Pd. Furthermore it is known that adsorbed 
hydrogen significantly reduces the inward relaxations at 
metal surfaces as it increases the bond coordination of 
the surface atoms, making them, to some extent, more 
bulk like. However, the possibility of hydrogen contami- 
nation was strongly rejected by later experimental papers 
(e.g. 10). 

Morrison et al. O] investigated an alternative possibil- 
ity ||ll| , namely that the presence of surface magnetism 
could increase the first interlayer spacing, i.e., reducing 
the large inward relaxation they had obtained in their 
non-magnetic calculation by "magnetic pressure" . In 
fact, bulk Rh is already close to fulfilling the Stoner cri- 
terion of ferromagnetism, and the narrower density of 
d-states at the surface might stabilize a magnetic state 
at the surface. Density-functional theory (DFT) together 
with the local-density-approximation (LDA) gives a non- 
magnetic ground state for Rh(OOl), but this might be 
due to the LDA. For example, for bulk iron, which is 
studied in greater detail, the LDA falsely puts the bcc 
magnetic ground state at a higher energy than the non- 



magnetic hep and fee states |l 
LDA problem Morrison et al. 



To get around this 



employed a pseudopo- 
tential which is based upon an atom in which all the 
electrons see a Hartree-Fock exchange potential arising 
from the core electrons and an LDA potential arising 
only from the valence electrons. Then in the surface 
calculations the valence exchange potential was taken 

1 /3 

proportional to "ra^aicnco- ^^ ^ consequence, they found 
that their Rh(OOl) surface is ferromagnetic. The mag- 
netic moment is M = 1.8 /is /surface atom, and the re- 
sulting magnetic pressure reduced the surface relaxation 
Adi2/do from —3.22 % (in the nonmagnetic equilibrium 
state) to —1.52 % in the magnetic ground state. Thus 
these authors concluded that surface ferromagnetism is 
the driving force giving rise to the small surface relax- 
ation deduced experimentally. Subsequently performed 
theoretical work, however, did not accept their approach 
and conclusions [||jlj] ; and also experimental studies pro- 
vided no convincing support |15[ . In their spin-polarized 
photoemission experiment Wu et al. |Iq] found only a 
weak indication of surface magnetism with a small mag- 
netic moment of about M — 0.2 /ig/surface atom. 

In this letter we present a new theoretical study which 
extends the previous work by considering the general- 
ized gradient approximation (GGA) |16|, and by taking 
zero-point effects and the thermal expansion as well as 
surface magnetism into account. Such a study is desir- 
able since all previous DFT calculations |jl|,||,||J^,^ were 
performed with the LDA which docs not describe the 
magnetic state reliably; furthermore, in all previous work 
zero-point and thermal vibrations were ignored, while the 
LEED data were taken at room temperature |p|,0J§] . We 
will show that the above noted discrepancy between the- 
oretical and measured results is mostly due to the un- 
justified neglect of vibrational contributions to the free 
energy. It is argued that the vibrational effects will typ- 
ically play a much bigger role than hitherto anticipated. 
Furthermore, we find that surface magnetism has a very 
small effect on the surface interlayer distance. 

We employ the full-potential LAPW method | |l7|,p^ to- 
gether with norm-conserving pseudopotentials W9 . The 
nonlinear core-valence exchange-correlation interaction 
is treated using the correct core-electron density as ob- 



tained in the atomic calculation ^^. The method gives 
an accurate and at the same time computationally ef- 
ficient description of the interatomic interactions, total 
energies, and stable or metastable geometries. Our GGA 
calculations are performed consistently by creating the 
pseudopotential from first-principles DFT-GGA calcula- 
tions. The Rh (001) surface is modeled by a periodic slab 
geometry consisting of nine layers of Rh and a vacuum 
thickness corresponding to five such layers. The geome- 
try is optimized by a damped molecular dynamics p8[ , 
allowing the top two layers on both sides of the slab to 
relax. The remaining atoms are kept at the bulk lattice 
sites. The parameters describing the LAPW basis set 
are: (K*^^)^ = 14 Ry and l^^^ = 8. For the k summa- 
tion we use 28 points of the irreducible part of the surface 
Brillouin zone. 

Since aU previous calculations |i|,|,|,|j| for Rh (001) 
were performed with the LDA, we also performed LDA 
calculations, which together with our GGA results al- 
low us to examine the effect of the GGA on the surface 
properties of Rh (001). Using DFT-LDA our bulk lattice 
constant is 3.81 A, which is in good agreement with pre- 
vious calculations P,p|. The experimental result, which 
unlike the quoted calculated value contains the influence 
of zero-point vibrations, is 3.79 A [^. Were the zero- 
point vibrations to be included in the theory, the calcu- 
lated lattice constant would increase by about 0.5 % p2[ . 

Using the GGA we find that the bulk lattice constant is 
expanded with respect to the LDA value by 2.2 %, givi ng 
it a value of 3.89 A. For hep Ru [|| and fee Pd Q a 
similar trend was found when comparing LDA and GGA 
lattice constants (see also Ref. [|5|). However, we find 
that the GGA affects the surface relaxation of Rh (001) 
only little (see Table I), although the cohesive energy, 
the surface energy, and the work function are affected 
noticeably compared to the LDA values. 

Table I summarizes the results for surface relaxations, 
work functions, and surface energies as obtained by dif- 
ferent calculations and experiments. With respect to 
the surface relaxation it is immediately evident that the 
LAPW calculations by Feibelman and Hamann ^ give 
an exceptionally large value. The present LDA calcula- 
tions, those of Cho and Kang [Q, and those of Methfes- 
sel et al. H, who did not relax the second layer, are in 
good agreement with each other. Also the result of the 
nonmagnetic study of Morrison et al. H (quoted above) 
agrees well with our value. As previously pointed out by 
Morrison et al. g, the too large top-layer relaxation in 
Feibelman and Hamann's calculations maybe attributed 
to the use of the poor k-point sampling [^ . 

The difference between our DFT-LDA results for the 
surface relaxation (A(ii2 /do = —3.0 %) and the previous 
LEED analyses IqJ^H is decreased significantly compared 
to the results of Ref. |],|9) (see Table I); the DFT-GGA 
calculations give a result [Adu/do — —2.8 %) which 
is even closer. We will now show that the physics of 
Rh(OOl) is much more interesting than previous studies 
had anticipated. At first we will address the influence of 



lattice vibrations of the Rh (001) surface and show that 
the restriction to the T — OK total energy falsely neglects 
some important physical aspects, which clearly affect the 
free energy and as a consequence the surface properties. 
Then we analyze the possibility of surface magnetism. 

It is well known that the zero-point vibrations give 
rise to a recognizable effect on the bulk lattice constant. 
Moruzzi et al. |2^ had systematically included this effect 
in their KKR studies of metals. Typically, however, this 
effect has been ignored. It is plausible that vibrational 
effects may be even larger at surfaces than in the bulk. In 
a correct treatment the equilibrium structure at a given 
temperature is determined by the minimum of the free 
energy. At not too high temperatures this differs from 
the total energy of the rigid lattice mainly by the con- 
tributions from atomic vibrations to the internal energy 
(including the zero-point vibrations) and the vibrational 
entropy. In the quasi-harmonic approximation the free 
energy for the surface is F{T) = Mind^^F{di2,T) with 



F{di2,T) = V{di2) + kBTY,{ 



+ln\ 1 — exp- 



huji{di2) 
2kBT 

Tiuj,{di2) 



knT 



■)} 



(1) 



where huji{di2) denotes the vibrational frequencies and 
the sum goes over all bands and k points. The first 
term in eq. (1) is the first interlayer potential and the 
second term is the vibrational energy and entropy. We 
note in passing that such a quasi-harmonic description 
had been used successfully in DFT calculations of the 
anomalous thermal expansion of covalent semiconduc- 
tors IpO]. For Ag, Cu, and Al surfaces eq. (1) has been 
recently evaluated by Narasimhan and Scheffler |31|. We 
note that the equilibrium distance di2 is shifted away 
from the minimum of V{di2) towards a larger interlayer 
spacing and that this shift is determined by the slope of 
the TiLJi{di2) but not their actual values. To a first ap- 
proximation this dependence of ui on di2 depends only 
weakly on the band index and k. We therefore replaced 
the sum in eq. (1) by three surface-phonon wave pack- 
ets. Only the top layer is moved, and deeper layers are 
kept fixed. Figure |l| provides our DFT-GGA result for 
the potential energy ^(^12); its curvature gives the fre- 
quencies for the perpendicular vibrational mode. For the 
the parallel vibrations we use two "modes" along [110] 
and [110], which are actually degenerate. The calculated 
phonon energies TiLOi of the in-plane and out of plane vi- 
brations are shown in Fig. 0. Our result for the tem- 
perature dependence Adi2{T)/do, considering the three 
above discussed phonon "modes" , is given by the full dots 
in Fig. y. It is obvious that thermal vibrations have in- 
deed a pronounced effect. They change the surface relax- 
ation from the value given by the minimum of the total 
energy, -2.8 %, to Adia/do = -1-4 % at 300 K. This 
result is now in excellent agreements with that of the 
room-temperature LEED analysis @] which determined 
a value of -1.16 ±1.6 %. 



It is interesting to note that the motion of the surface 
layer parallel to the substrate yields the most important 
contributions (compare Ref. 1^^). If we would neglect 
the contributions of the parallel motion and use only the 
perpendicular vibration the resulting top-layer relaxation 
would be much smaller. This result, displayed by the 
open dots in Fig. 0, reveals that the anharmonicity of 
the interlayer potential of Rh (001) does not have a very 
pronounced influence on the top-layer relaxation. 

Our DFT-GGA calculations predict that the ground 
state of Rh(OOl) is nonmagnetic. This result remains 
even if we intentionally increase di2 to the unrelaxed 
geometry, thus offering a bigger volume per surface 
atom which typically helps to stabilize a magnetic state. 
Despite this apparently clear result of a non-magnetic 
ground state, we asked how far away in energy the fer- 
romagnetic state actually might be. For this purpose 
we performed spin-polarized calculations employing the 
fixed-spin- moment method [^ . Figure Q shows the total 
energy versus magnetic moment for a given relaxed sur- 
face of Adi2/do = —2.4 %. We find that the total energy 
monotonically increases with increasing magnetic mo- 
ment. This behavior is similar to that of a previous fixed- 
spin-moment study of Cho and Kang [0, who used the 
LDA. The present DFT-GGA result for the energy dif- 
ference AE between the nonmagnetic and ferromagnetic 
states is however reduced significantly compared to the 
previous LDA one jj], and Fig. Q shows that AE remains 
almost constant until the magnetic moment reaches a 
value of 0.5 /i^/surface atom [ p3| . In the fixed-spin- 
moment method p3| , spin-up and spin-down eigenvalues 
are calculated for different Fermi energies. For a mag- 
netic moment of M — 0.5 /is/surface atom we find that 
the difference between the two Fermi energies is only 25.9 
meV; the total-energy difference at M = 0.5 /is/surface 
atom is only 1 meV. In other words, our calculations 
show that the ferromagnetic state is practically degen- 
erate with the nonmagnetic one, and we expect that a 
weak ferromagnetic state will occur on Rh (001) possibly 
stabilized by surface imperfections. This result is consis- 
tent with the room temperature spin-polarized photoe- 
mission experiments by Wu et al. |15|, who observed a 
rather weak ferromagnetism with the surface magnetic 
moment of about 0.2 /zs/surface atom. To some extent 
our results support the motivation behind the study of 
Morrison et al. 0, although their treatment predicted a 
rather strong and stable ferromagnetic state. In contrast 
to them we find that the magnetic state is very close to 
the critical point, that the magnetic moment should be 
very small, and thus the magnetism has practically no 
effect on the surface relaxation, or vice versa. 
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TABLE L Surface relaxations Adi2/do and Ad23/do (do 
is the bulk interlayer spacing), work functions (eV), and 
surface energies a (eV/atom) for Rh (001) as obtained by dif- 
ferent calculations and experiments. 



Adi2/do 



Ad23/do 



LDA |1 
LDA [2 
LDA |i 
LDA I 
this - LDA 
this - GGA 



-5.1% 
-3.5% 
-3.8% 
-3.8% 
-3.0% 
-2.8% 



-0.5% 



-HO.7% 
-0.2% 
-0.1% 



5.49 
5.25 



5.26 
4.92 



experiments -^0.50 ± 1.0% [7| ± 1.5% 7| 4.65 [21 



experiments -1.16 ±1.6% [8| ± 1.6% 



8| 4.98 [27 



1.12 
1.27 
1.29 
1.44 
1.29 
1.04 
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FIG. 1. Total energy per surface atom as a function of the 
top-layer relaxation for Rh(OOl). The minimum of the fitted 
curve is set to be the energy zero. 
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FIG. 2. Phonon energies tiuji of the in-plane (solid line) and 
out-of-plane (dashed line) modes of Rh (001) as a function of 
the top-layer relaxation. 
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FIG. 3. Top-layer relaxation of Rh (001) as a function of 
temperature. Full dots represent results obtained using eq. 
(1) with the results of Figs, hi and bl Open dots show results 
obtained if the parallel vibrations are neglected. 
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FIG. 4. Total energy per surface atom as a function of the 
magnetic moment per surface atom for a surface relaxed by 
Adi2/do = — 2.4 %. The nonmagnetic state defines the energy 
zero. 



